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ABSTRACT 

We study the variations in the Cyclotron Resonant Scattering Feature (CRSF) during 
2011 outburst of the high mass X-ray binary 4U 0115+63 using observations performed 
with Suzaku, RXTE, Swift and INTEGRAL satellites. The wide-band spectral data 
with low energy coverage allowed us to characterize the broadband continuum and 
detect the CRSFs. We find that the broadband continuum is adequately described 
by a combination of a low temperature {kT ^ 0.8 keV) blackbody and a power-law 
with high energy cutoff (ifcut "^5.4 keV) without the need for a broad Gaussian at 
^ 10 keV as used in some earlier studies. Though winds from the companion can 
affect the emission from the neutron star at low energies (< 3 keV), the blackbody 
component shows a significant presence in our continuum model. We report evidence 
for the possible presence of two independent sets of CRSFs with fundamentals at 
^11 keV and 15 keV. These two sets of CRSFs could arise from spatially distinct 
emitting regions. We also find evidence for variations in the line equivalent widths, 
with the 11 keV CRSF weakening and the 15 keV line strengthening with decreasing 
luminosity. Finally, we propose that the reason for the earlier observed anti-correlation 
of line energy with luminosity could be due to modelling of these two independent line 
sets ('^11 keV and 15 keV) as a single CRSF. 
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1 INTRODUCTION 


4U 0115+63 is a high mass X-ray binary s ystem, first dis¬ 
covered in the UHURU sat ellite’s sky survey dCiacconi et al.l 
I 1972 I : [Porman et al.lll97^ 1. with more than 1 5 subsequent 
outbursts recorded till date (iBoldin et al.llioi^ l. The system 
consists of a pulsating neutron star with spin period ~ 3.61s 
( Cominskv et al.l 1978l l and a B0.2Ve main sequence star 
( Johns et al.lll9780. wi th an orbital period of ~ 24.3 days 


i Rappaport et al]ll978l '). The distan ce to this binary system 
has b een estimated to be ~ 7 kpc (iNeeueruela fc Okazakil 
l200lh . The source exhibits luminous Type II X-ray out¬ 
burst during which multiple cyclotron resonance scatter¬ 
ing features (CRSF) have been observed in the X-ray spec- 
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trum, with 5 detec ted harmonics llSantangelo et ahl Il999l : 
iFerrigno et al.l[2009l l. CRSF are caused by scattering of X- 
ray photons from electrons in the accreting plasma chan¬ 
neled by the NS magnetic field. The energy at which these 
lines occur is given as Ecyc = A.GB 12 x (1 + 2 )“^ keV 
dCoburn et ahlfioo^ . Here B ^2 is the local magnetic field 
(in units of 10^^ Gauss) and 2 is the gravitational redshift 
in line energy. Thus cyclotron lines give us a direct probe of 
the local magnetic field near the scattering regions. 

Cyclotron line parameters of many sources are found 
to vary with the phase of rotation (see iHeindl et al.l |2004| . 
for a review), and the varying luminosity of the out¬ 
burst dBecker et al.l l2012l l. However, the variation in en¬ 
ergy of the fundamental CRSF of 4U 0115+63 with lumi¬ 
nosity has been the source of some debate. iNakaiima et al.l 
d200^ : [Tsvgankov et al.l d2007ll : lLi et al.l d2012ll find an anti- 
correlation b e tween the line energy and luminosity, whereas 
iMfiller et al.l d2013h find this anti-correlation to be an ar- 
tifact of the continuum spectral modelling. iBoldin et al.l 
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Figure 1. Variations in the count rate and luminosity during 
the 2011 outburst of 4U 0115+63. The continuous line repre¬ 
sents the observed count rate measured with the MAXI sky mon¬ 
itor. The symbols represent the pointed observations and the 
corresponding 3 — 50 keV luminosities derived from our spec¬ 
tral analysis. The two simultaneous wide-band observations are 
annotated separately. The horizontal line shows the luminosity 
level(4 X 10^^e rgss~^) around which the previous observations 
jLi et al have found a sharp change in fundamental CRSF 

energy. 


l|2013h have pointed out the cause of this reported dichotomy 
to be due to the use of a broad Gaussian like emission feature 
to model the continuum in some of the works. 

In this paper, we study multiple observations during 
the 2011 outburst of 4U 0115+63 performed with Suzaku, 
RXTE, Swift and INTEGRAL satellites, providing us with 
a wide-band coverage from 0.5 keV to ~ 60 keV with high 
signal-to-noise. As shall be seen in subsequent sections, this 
availability of wide bandwidth data is very important to 
correctly model both the continuum and the CRSF in the 
source spectrum. From the results of our spectral analysis, 
we find evidence for two sets of cyclotron lines whose param¬ 
eters vary with source luminosity. In the following sections, 
we describe the observations and reduction of data that we 
used in spectral analysis of these data-sets in ^ and the 
inferences and possible implications in ^ 


2 OBSERVATIONS AND DATA REDUCTION 

In all, we analysed data obtained by different X-ray obser¬ 
vatories over the 15 days of the 2011 outburst. Table [T] lists 
all the RXTE and Suzaku observations which were available 
in the archive and a Swift and an INTEGRAL observation, 
made when the source was near its peak luminosity. Figure[T] 
shows the variations in the count rates as measured with 
MAXI sky monitor during the 2011 outburst. The pointed 
observations performed at different luminosity levels are also 
marked in this figure. The reduction of data-sets from each 
of these satellites is explained in the following sub-sections. 


Table 1. A summary of pointed observations during the 2011 
outburst of 4U 0115+63. 


MJD 

Instrument 

Obsid 

MAXI“ 

Exposure^ 

55730.06 

RXTE 

96032-01-01-00 

0.788 

6.53 

55736.34 

Swift 

00031172010 

0.925 

6.69 

55736.34 

INTEGRAL 

106100650010 

0.925 

2.13 

55738.87 

RXTE 

96032-01-02-00 

0.993 

4.52 

55739.51 

RXTE 

96032-01-02-01 

0.745 

3.69 

55743.91 

RXTE 

96032-01-03-00 

0.682 

5.02 

55746.90 

RXTE 

96032-01-03-02 

0.526 

9.04 

55747.53 

Suzaku 

406048010 

0.526 

24.27 

55750.82 

RXTE 

96032-01-04-00 

0.352 

16.91 

55750.82 

Suzaku 

406049010 

0.352 

81.68 

55751.32 

RXTE 

96032-01-04-02 

0.388 

2.08 

55753.00 

RXTE 

96032-01-04-03 

0.248 

1.85 

55753.42 

RXTE 

96032-01-04-04 

0.238 

0.73 


“MAXI rate (in counts cm ^ s is of nearest observation. 
^Exposure time is in kilo-seconds. 


2.1 RXTE Observations 

Among the RXTE observations, we used all but two of them. 
Obs-Id 96032-01-04-03 and 96032-01-04-04 were made when 
source flux was very low, and with short exposure times 
(1856s and 736s respectively). Due to poor count statistics, 
we were unable to constrain the CRSF parameters for these 
observations, making them unsuitable for the present work. 
RXTE spectra were obtained from the raw data files us¬ 
ing ETOOLS from HEASoft v 6.15.1. We used data from 
both the HEXTE and the PGA detectors. PCU2 Science 
Array data were used for generating PGA spectrum, and 
Cluster A Science Array data were used for HEXTE spec¬ 
trum. HEXTE background was obtained from Cluster B, 
and dead-time corrections were applied to both the source 
and background HEXTE spectra. For PGA, we found the 
dead-time to be a maximum of 5%, which lowered its flux by 
about 5%. We did separately correct for this. Data grouping 
and usage of systematic errors for the RXTE spectra are 
shown in Tabled In all, the RXTE provided usable data 
covering the 3 keV to 50 keV band. 


2.2 Suzaku Observations 

Suzaku had two long duration observations made during 
the course of the 2011 outburst. For reduction of Suzaku 
data, we used ETOOLS from HEASoft 6.15.1 with CALDB 
updated till June 2014. We used data from the XIS ( 0.6 
keV - 10 keV ) and PIN ( 15 keV - 60 keV ). The 
GSO data had very low SNR and hence was not used. Al¬ 
though the PIN does collect data from 12 keV onwards, 
we discarded the data from 12 keV to 15 keV due to 
high uncertainty in PIN backgrouncG. XIS data above 10 
keV and PIN data above 60 keV were too noisy to be of 
use for our analysis. As a result we were not able to effec¬ 
tively cover the energy range (10 keV - 15 keV ), which is 
important to model the fundamental cyclotron line in this 
system. To overcome this, we used simultaneously taken 


see Suzaku Data Reduction ABC guide (version 4.0) 
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RXTE observations in conjunction with our Suzaku data¬ 
sets. Since obs-id 4060490 LO (Suzaku) and Obs-Id 96032-01- 
04-00 (RXTE) had some overlap, we used the portion where 
overlap existed. This was done by creating a GTI which cov¬ 
ered the common interval between the RXTE and Suzaku 
GTI hies. However Obs-Id 406048010 (Suzaku) had no over¬ 
lap with any other RXTE observation. Thus, we did not use 
this data-set. 

The XIS data were taken in the 1/4 windowed mode 
(of size 256 pixels) at normal clocking speeds to reduce the 
effect of pile-up. We found signihcant pile-up upto a max¬ 
imum level of 17% in the central regions and corrected 
for it using the recipe of John Davifl This was achieved 
by rejecting a central circular region of size ~ 25" from the 
XIS image for which the computed Pile-up percentage was 
greater than 6%. After extraction of data from individual 
XIS chips, the XIS 0 and 3 data were combined using the 
addascaspec tool. This resulted in two sets of spectral hies; 
one from the back illuminated (BI) GCD and the other from 
the combined front illuminated (FI) CCDs. XIS data below 
0.8 keV were rejected owing to discrepancies between the 
FI and BI spectra at a 3cr level. A mismatch at a ~ 2 (t 
level between the two spectra around the Si K-edge was also 
noticed, but we did not discard this part of the spectrum as 
it was important for htting our continuum model. This led 
to poorer values of the statistic that we use for htting 
the spectral data. 

The PIN data were extracted using the tool 
hxdpinxbpi. In the joint analysis of Suzaku and RXTE, 
we required a cross-normalisation of ~ 1.57, between the 
Suzaku XIS and PIN, much larger than the recommended 
value of 1.16 given in the Suzaku Data Reduction ABC 
guide (version 4.0). Therefore, we investigated this issue in 
some detail. To verify the PIN data, we extracted individ¬ 
ual 64 PIN count spectrum as explained in the Suzaku note 
on estimating P/A noise 0, and found no discrepancies. We 
also extracted night-earth data for the PIN observations to 
compare and see if the estimated backgrounds were correct. 
Fig. [^compares the night-earth and the “tuned” background 
spectral data. The two background spectra are similar in 
shape and in hux. 

We further examined the validity of the background 
subtraction by comparing the net (source minus back¬ 
ground) and background light curves. The lack of correlation 
amongst the two indicates that the PIN background was es¬ 
timated correctly. Finally we checked if the XIS spectra were 
properly area corrected. The XIS data were taken in the 1/4 
window mode with normal clocking, thereby making a rect¬ 
angular source footprint. By choosing both rectangular and 
circular source regions for extraction we found no difference 
in the resulting spectrum. The difference in the fit-quality 
as a result of using the recommended and the best-fit cross¬ 
normalisation factor, 1.16 and 1.57 respectively, is shown in 
Figure [21 As seen in the figure, the discrepancy in the ra¬ 
tio plot for Suzaku PIN (seen as a set offset points in the 
ratio plot of the left panel) vanishes when we use the best 
fit cross normalisation factor. We have thus used a cross- 


^ http://space.mit.edu/CXC/software/suzaku/ 

® http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/pinnoise.html 


normalisation factor of 1.57 for all our subsequent models 
and fits (also see Sect. 3). 

2.3 Swift and INTEGRAL observation 

We selected this particular observation from multiple Swift 
observations of the source, because it was simultaneously 
taken with an INTEGRAL observation (thereby enabling 
wide-band spectral coverage) and it was near the point hav¬ 
ing highest MAXI counts. This observation gave us a wide¬ 
band data-set near the peak luminosity with similar energy 
coverage to the Suzaku data-set at much lower luminosity 
levels. However, as pointed out in section [S] the low effec¬ 
tive area of INTEGRAL JEM-X resulted in lower signal to 
noise in the energy band from 11 keV to 17 keV , than for 
the Suzaku data-set. 

Swift data reduction was carried out using ETOOLS 
from HEASoft 6.15.1 with CALDB updated till June 2014. 
INTEGRAL analysis was done using OSA 10.0 with its cal¬ 
ibration files also updated till June 2014. We used the Swift 
XRT data from 0.5 keV to 9.7 keV for our spectral anal¬ 
ysis. These data were taken in the windowed timing (WT) 
mode of XRT. The maximum observed count rate was about 
40 cou nts s~^. This ensured that the XRT data was not piled 
up fsee lRomano et ^l2006l l. The XRT data from 0.4 keV to 
1 keV did show an excess due to uncertainty in the response 
modelling of XRT WT mode, but upon using the position 
dependent WT response files, this apparent excess was re- 
movec0. Hence we used these position dependant response 
files for all our analysis. 

We used two instruments - the JEM-X and the IBIS 
from the INTEGRAL’S suite. The JEM-X 2 data (from 
JEM-X) and the ISGRI data (from IBIS) were extracted 
with higher spectral binning than used in the standard 
pipeline. The JEM-X data were extracted from 3 keV to 
30 keV and the IBIS data from 15 keV to 100 keV , as per 
the instructions in the IBIS Analysis User Manual (Issue 
10.0) . The latest calibration files were used for generating 
the spectra and the rebinned rmfs. After extraction, the data 
were grouped and systematic errors added as given in Ta¬ 
ble [5] We used spin phase averaged spectra in all the obser¬ 
vations. Foimnodelling the spectra, we used XSPEC vl2.8. Ig 
(lArnaudlll996l ~l. 


3 SPECTRAL ANALYSIS 

We began our spectral analysis with the wide-band Suzaku 
dataset (obs-id 406049010) taken alongwith RATE dataset 
(obs-id 96032-01-04-00). For other datasets, we used the re¬ 
sults from this analysis as a template because this data-set 
had both wide-band spectral coverage and high signal to 
noise. As noted previously, when we held the PIN-XIS cross¬ 
normalisation fixed at 1.16 we found a significant offset of 
only the PIN spectrum from the other instrument spectrum 
in the ratio plots (refer Fig. [3] for the best fit cutoffpl based 
model with CRSFs). When we let the normalisation param¬ 
eter free, we found that it gave a best fit value of 1.57 which 
we adopted for all our subsequent analysis. 

http://www.swift.ac.uk/analysis/xrt/rmfs.php 






4 Iyer et al. 



Correlation ourve 



-3.9E+04 -1.3E + 04 1.3E + 04 3.9E+04 

Shift {§ Bins) 


Figure 2. A comparison of the estimated tuned PIN background and the night-earth PIN data. The left panel compares the spectral 
data extracted from reconstructed PIN background and the Night Earth data. See online edition for a color print of this figure. Right 
panel show the cross-correlation between the net (source - background) and background lightcurves. See text for details. 
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Figure 3. The spectral data {XIS data in red and black, PIN data in dark blue, PCA data in green and HEXTE in light blue color) 
and the ratio plots for best-fitting cutoffpl model. These plots were made with the A/S'-to-PJA cross-normalisation fixed at 1.16 (left 
panel) and with the best-fit value of 1.57 (right panel). The offset is seen in the ratio plot in left panel when using a cross-normalisation 
value of 1.16. See beginning of section 3 for comments on this figure and the online edition for a color print of this figure. 


Table 2. Grouping scheme and systematic errors used for differ¬ 
ent spectral datasets. 


Spectral data 

Grouping 

scheme 

Systematic 

errors 

RXTE PCU 

None 

1% 

RXTE HEXTE 

Minimum 60 counts/bin 

1% 

Suzaku XIS 

Minimum 60 counts/bin 
and oversample by 3 

None 

Suzaku PIN 

Minimum 40 counts/bin 

None 

Swift XRT 

Minimum 60 counts/bin 

None 

INTEGRAL JEM-X “ 

31 bins in 3-30 keV 

3% 

INTEGRAL IBIS “ 

27 bins in 15-100 keV 

1% 


“ Coded mask datasets have Gaussian and not Poisson statistics 
and therefore do not need rebinning. 


effects of local absorption and emission features from the 
plasma wi nd of the co mpanion star are known to affect the 
spectrum (|Suc However, we did not specifically ac¬ 

count for this in the spectrum. Complete modelling of the 
wind spectrum is preferably done with higher spectral res¬ 
olution data from gratings. As mentioned in Sec. 2.2, the 
Suzaku data-set had calibration uncertainties around the 
detector’s Si K-edge (from ~ 1.6 keV to 2.5 keV ). Thus 
modelling the wind effects was difficult with this data-set. As 
these effects are restricted to lower energies, and as we do not 
consider data below 0.8 keV in our analysis, we do not expect 
this to change our continuum model estimates significantly. 
The interstellar absorption was mo delled by an updated ver¬ 
sion of tbabsQ using abu ndances oflW ilms et ^ (l2000l f and 
cross-sections as given in IVerner &: Yakovlev ( lOOsI l. 

This model fit well if we took data above 3 keV , but 
showed a significant excess when we included the data from 


To model this wide-band spectrum, we started off by 
using the simple cutoffpl based continuum modified by inter¬ 
stellar absorption and CRSFs. At low energies ( < 2 keV), ® http://pulsar.sternwarte.uni-erlangen.de/wilms/research/tbabs 
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0.8 keV to 3 keV . This is clearly seen in the top left panel of 
FigurelH where we obtained /dof — 2018.60/474. On us¬ 
ing a blackbody (bbody) to account for this excess, we fonnd 
that the spectral data were well fitted to give a significant 
improvement in fit /dof = 731/472. The temperature and 
radius of the blackbody so obtained are listed in Table|31 The 
radius was computed from the normalisation of the black¬ 
body model bbody, which depends on the luminosity and 
distance to source. The distance was taken as 7 kpc and 
luminosity was taken as Lbb = cr. T'^.Area. The blackbody 
radius is clearly a factor of 10 larger than expected values 
for a hotspot at the base of the accreting mound. However, 
given that we might be integrating emission area from the 
base of the column and that effects of comptonization of this 
blackbody due to the accreting plasma are not considered, 
it is likely that these val ues are over-est i mates . We did try 
the XSPEC model given bv iFarinelli et all (l2012lf to see if this 
was true and found values of blackbody radius to drop to 
about 6 kms (see Appendix|^for details). The large number 
of free parameters in this model, though, make it difficult 
to use while testing the CRSF parameters across multiple 
data-sets. 

We also tried using comptt model to account for t his low 
energy excess, as indeed was tried in iFerrigno et al] (120091 1. 
This gave us an improvement in fit with the reducing 
from 688/469 to 583/466. This large improvement strongly 
indicates that the blackbody seed photons are indeed comp- 
tonized. The spectral fits for comptt were obtained though, 
for a electron plasma temperature (kTe) of 1.2 keV , a seed 
photon temperature of 0.2 keV and plasma optical depth(r) 
of 35. The comptt code itself is meant for use for plasma tem¬ 
perature greater than 2 keV and at such high r, the shape of 
the spectrum is very similar to the Wien tail of a blackbody. 
Using the comptt model for this Wien tail alone, also lead 
to an increase in the value of the absorption column from 
Nh = 1.3 X 10^^ cm“^ to Nh = 1.7 x 10^^ cm“^. When we 
tried the same model in the Swift / INTEGRAL spectrum, 
we found a very small reduction in the ht statistic x^ from 
146/166 to 145/163. Furthermore, none of the RXTE obser¬ 
vations gave an improvement in fit x^ when using the comptt 
model as compared to the bbody model. Finally, we found 
very little differences in the CRSF energies, while using the 


Table 3. Blackbody radius and area for different observations. 
Note that RXTE observations have higher errors. 


Instrument 


Day (MJD) kT ( keV Radius (kms;[^ 


RXTE 

Swift / INTEGRAL 

RXTE 

RXTE 

RXTE 

RXTE 

Suzaku / RXTE 
RXTE 


55730.06 

1 

-*-•^^+0.19 

13 

55736.34 

0 

^•^‘^+0.41 

1 7 3Q“‘^-^^ 

± t •'3i?+9.22 

55738.87 

U.y4_^0.65 

1 1 QS — 6.65 

-L < •'30^5^5 3g 

55739.51 

i.uy_l_Q 0 

14 69“^'® ^ 

iq:.Dy_|_Q Q _ 

55743.91 

o-79;S:2? 


55747.9ljSl 

o-78;n 

18.85;°;° 

55750.82 

<'3+0.01 

13.95;°:« 

55751.32 

u.y^_l_0 38 

10 54“^'^^ 

1U.04_^5 37 


“ All errors in this table are la deviations 
^ calculated as 

Could not put error constraints on blackbody temperature 


Table 5. Continuum parameters for the two wide-band observa¬ 
tions 


Parameter 

Swift / INTEGRAL 

Suzaku / RXTE 

Nh( 1022 cm-2) 

1.3 

1.3 

r 

^•-‘-’^+0.48 

. ^c-O.OS 

-L-^'3_,_q Q7 

Ecut keV 


5 34“^-^® 

'3-0^_,_q 21 

keV 

'-'•^^+0.41 

0 73“0 0l 
'J- '-J+o.oi 

bbnorm (x 10 

-*-^•^+1.6 

o y —0.2 
'^•'+0.2 

Fe line eq.width ( eV) 

-^-^•'+20.5 

11 7-3.0 


improvement in /dof from 756/473 to 731/472. Upon 
testing this with the simftest script of XSPEC, we found it 
to be significant to a level greater than 3 (t for 480 itera¬ 
tions. This test was done using XIS data only, to isolate any 
effects of the continuum beyond 10 keV . 

We model the CRSFs us ing a symmetric Lor entzian like 
absorption feature cyclabs (iMihara et al.lll990l l. with each 
line energy being independently determined by the fit. This 
was done to account for anharmonic line ratios which can oc¬ 
cur due to effects of viewing geometry and asymmetric emis- 


simple blackbody model as compared to the comptonized 

sion 

patterns, as has been reported in other HMXB sources 

blackbody. Since the data that we analysed could not help 

(iKrevkenbohm et al.ll2002l: Pottschmidt et al.l 2005h and in 

us pin down the nature of the low energy continuum, we 

theoretical simulations dNishimural l2005l: ISchonherr et al. 

used the simple blackbody for all subsequent analysis. How- 

2007 

: iMukheriee & Bhattacharval I2OI2I: iMukheriee et al. 


ever, we do note that probing this component further may 
be important to understand the overall continuum, which in 
turn may affect the cyclotron line results. 

As stated in ISuc ii^ (l201l[l and iMfiller et al.l ll2013ll , we 
found the column density (Nh) to be strongly correlated 
with the power-law index (F) and to be varying across ob¬ 
servations. To prevent any unwanted effects due to changing 
Nh on our fits, we fixed the column density for all observa¬ 
tions to the best fit value of Nh ~ 1.3 x 10^^ cm“^ obtained 
from the wide-band Suzaku observation. We also used a nar¬ 
row emission line with its centroid fixed at 6.4 keV and width 
fixed at lO”"^ keV to model the Fe K-a fluorescence emis¬ 
sion which has been observed in many other HMXB NS bi- 
nary systems and is also expected to be seen in 4U 0115-1-63 
llTorreion et ahll^OlOl : IMfiller et al.ll2013l l. For the wide-band 
Suzaku observation, this narrow Fe-line component gave an 


201311 . 


Modelling the Suzaku spectrum, we found cyclotron 
lines to be at ~ 14 keV , ~ 21 keV and ~ 31 keV . How¬ 
ever, we additionally noted an absorption feature at 11 keV 
in the fit residuals. Modelling this using a Lorentzian like 
absorption feature gave a significant improvement in the fit 
as seen in Fig. [4] with the /dof changing from 731/472 
to 688/469. We were not able to ht for any CRSF fea¬ 
tures above 35 keV in any of the datasets because of large 
noise at the higher energies. And the CRSF at ~ 33 keV 
was itself constrained by keeping its FW HM hxed at 4 keV 
dFerrigno et al.l[2009l : iMuller et cd1l2013ll . The 4 CRSFs ob¬ 
tained in this observation (as listed in Table [4]) are at about 
11 keV , 15 keV , 20 keV and 33 keV . 

We tried next to model the high luminosity wide-band 
data of Swift and INTEGRAL simultaneous observation. We 
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Figure 4. Results of our joint spectral fits to the Suzaku XIS (black and red), PIN (blue) and RXTE PCA (green) and HEXTE (light 
blue) data using a cut-off power-law, CRSFs and a narrow iron Ka line at 6.4 keV. Top left: Unfolded spectral data, the model fitted 
above 3 keV and extended to lower energies in the upper panel and the deviations of the data from the model in the lower panel showing 
an excess at lower energies. Top right: Similar to top left but the model includes a low temperature blackbody in addition to the cutoffpl 
and 3 CRSFs. The residuals show possible presence of a 4th cyclotron line at ~ 11 keV. Bottom left: Similar to top right but the model 
includes a 4th cyclotron line at 11 keV . Bottom right: The 4 cyclotron lines as seen in the ratio plots obtained by setting the optical 
depth of the lines to zero after obtaining the best fit. See online edition for a color print of this figure. 


were able to fit the spectrum from 0.5 keV upto 60 keV us¬ 
ing only the cutoffpl continuum with absorption fixed at 
Nh = 1.3 X 10^^ cm~^ and CRSFs ( 11 keV , 20 keV and 

33 keV ) to get a /dof — 172.8/168. On adding the black- 
body component, this improved to fdof = 146.1/166. We 
tested the significance of this using the Irt script of XSPEC. 
This gave a significance of greater than 3cr for the blackbody 
model for 408 iteration runs. When taken alongwith the case 
of Suzaku / RXTE data-set, we see that the presence of a 
soft X-ray component is justified. As stated above, model¬ 
ing this component is important. For the sake of consistency, 
we use the blackbody model, as using the comptt model for 
the soft X-ray component did not give us any improvement 
over the blackbody based model for this observation. While 
using this blackbody and cutoffpl continuum, we could not 
fit this data-set with the 5th line as earlier and the fit was 
consistent with only 3 CRSFs. A summary of the contin¬ 
uum parameters obtained for the two wide-band data-sets 
is given in Table [5] Errors, unless other-wise mentioned are 
quoted for a level of 90% confidence in all cases. 

We finally analysed the set of standalone RXTE obser¬ 
vations. Due to the lack of coverage of the low energy bands 
in these data-sets, it was difficult to constrain the blackbody 
parameters. Hence we fit the blackbody by starting from a 


guess value based on the fit parameters obtained from the 
wide-band data-set and let the fit routine converge to give 
the values as quoted in Table [3] The blackbody tempera¬ 
ture was frozen to this best fit value while computing errors 
on other fit parameters. Excepting one, all the RXTE data¬ 
sets gave an improvement in the fit statistic when using 4 
CRSFs as compared to the fit with 3 CRSFs (see Table |4]). 
The exception was the first data-set, taken during the rising 
phase, which needed only two CRSFs to describe the spec¬ 
trum. We consider this case separately in the Sec. 4. Given 
this improvement when using 4 CRSFs, we decided to test 
its statistical significance for describing the spectral data. 


3.1 Statistical significance of the 4*^ cyclotron line 

Computing detection significance of CRSFs must be handled 
differently from the standard techniques usually employed 
for emission features. F irst, following the method outlined in 
lOrlandini et al.l (l2012l ) we construct the F-statistic directly 
as a ratio of the normalised fit dOrlandini et al.l l2012l : 
iBevington fc Robinson|[r992l ). as against using the standard 
F-test implemented in XSPEC, which constructs the F-stat as 
a ratio of change in normalised fo the original y^. 
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Table 4. Best-fit parameters for the CRSFs. 


Day (MJD) 

Ecyc ( keV ) 

FWHM ( keV ) 

r 

X^/dof 

Eq. width ( keV ) 

Ax^ / Ado^ 

F-stat % 

55730.06 

XU.OO^Q gg 

^•'^^+0.56 

-‘-•^'J+0.16 

57 

/ 

73 

11.98 ± 0.72 

0/0 

0 

(RXTE) 

37.29;}:J| 

4.0 +?,0 

^•^^+0.14 




2.38^ 



55736.34 

10 31 “4.58 

5 Qg“3-63 
0.yo_,_4 58 

0 97“0 ®2 

'^•^'+0.31 

146 

/ 

166 

6.21 ± 4.23 

0/0 

0 

{Swift / INTEGRAL) 

2o.oo;®:2i 

^•^^+5.24 

-*-•^^+0.31 




11.98 ± 4.50 




35 75 — 2.19 
Ci0./0^4 17 

4 

^•^+0.0 

^•^^+0.38 




2.3^^ 



55738.87 


4 33“1-47 
^•'^'^+1.22 

0 62“^’^ 
^•'^^+0.25 

47 

/ 

67 

3.35 ± 1.18 

14.91 / 3 

83 

{RXTE) 

15.82;0:«« 

-J .-1.06 

-*-•-*-+2.97 

n 1 —0.06 
'^•-*-+0.31 




0.17 ± 0.19 




19.86;?;®? 

Q Qo — 1.62 

^■°^+1.91 

-^•^-^+0.38 




10.48 ± 1.59 




35 . 53 ;?:?® 

4 0~^ 

^•^+0.0 

0 35“^-^^ 

'^•'50+0.16 




1.9i^ 



55739.51 

11 66“^'^^ 

11.00^0.39 

5 1 9 “^-^ 
^•-*-^+1.98 

n 74-0-44 
'^+0.2 

43 

/ 

67 

4.49 ± 4.56 

5.81 / 3 

64 

{RXTE) 

16.08;?:?? 

^•^^+4.53 

0 14”^’^ 
^•^^+0.98 




0.42 ± 1.40 




20 . 42 ;?:?? 

^■^^+1.25 

1 01 -0.22 

-‘-•'^-*-+0.13 




10.58 ± 3.59 





4 

^•^+0.0 

0.38;?:?? 




2.0iEl 



55743.91 

ll.zo^l Q3 

4 13 — 2.37 
^•-^'^+5.04 

0'45;?:?? 

31 

/ 

67 

2.50 ± 2.94 

10.17 / 3 

84 

{RXTE) 

15 . 2 ;?:?? 

0 89 — 

^•°^+3.08 

0-38;?:?? 




1.44 ± 2.35 




19.56;?:?? 

7 57“1’^3 
'•^'+1.39 

1 IJ-0A2 
-^•-^'+0.11 




8.86 ± 2.09 




32 . 77 ;?:?? 

4 0~^ 

^•^+0.0 

^•'^+0.15 




2.9iEl 



55747.90 

ll.ZU^l Q4 

0 50 — 2.01 

Q 00 -U.iy 

52 

/ 

67 

1.63 ± 1.11 

14.99 / 3 

81 

{RXTE) 

14 . 93 ;?:?? 

2 76“^'^ 

'’^+2.41 

0 35 — ^-^® 
'^•'50+0.51 




1.33 ± 1.40 




18 . 01 ;?:?? 

^•^^+1.67 

-| 00 — 0.32 
-*-•'^‘^+0.1 




11.67 ± 1.34 




33 7q-l-96 

00 .1 y+2.52 

4 0~® 

^•^+0.0 

^•^^+0.15 




1.4^ 



55750.82 

10 77“0-27 
'+0.32 

'^•^^+0.67 

0 12“0U3 

'^•4^4-0.04 

688 

/ 

469 

0.17 ± 0.10 

42.68 / 3 

74 

{RXTE / Suzaku) 

14 39“0-33 

i^.oy^Q 35 

4 19 —o.y:^ 

^•-Ly+ 0.97 

0 69~o-^^ 

u.uy_l_0 22 




3.48 ± 1.60 




ly.oo^l 4g 

Q 00 — 1 • 52 

^■'^'^+1.56 

-^•-^'^+0.14 




9.52 ± 2.21 




31 . 35 ;?:?? 

4 0~® 

^•^+0.0 

O.OD_^0.16 




2.8iil 



55751.32 

10.76;?,o 

2 l“2-l 
^•-^+7.33 

^•-^^4-0.59 

57 

/ 

67 

0.37 ± 0.59 

4.32 / 3 

55 

{RXTE) 

I4.5i;?:?? 

4 67“2-91 
^•^'+3.08 

0 73“0-28 
'*^4-0.61 




4.06 ± 1.07 




19.55;?:?" 

0 ij — 1.62 

°-'+13.52 

1 34~0-62 




10.90 ± 1.84 




31.98;?:?? 

4 0~^ 

^•^+0.0 

1 05~^’^ 
-*-•^^4-0.42 




4.4/0 




“ Change / reduction in on addition of extra CRSF / change in number of dof - denotes the confidence that a 4 CRSF model is 
better than a 3 CRSF model 
^ Cannot quote error on this, as width is frozen 


We try to compute the significance of the 4th CRSF 
at 11 keV , which we used as stated previously to improve 
the fit. For the RXTE/Suzaku joint spectral fit, the above 
method gave us an F-statistic value of 1.06 with a probabil¬ 
ity value for occurrence due to random noise (PCI) equal to 
0.26. We also noted that ignoring data from 1.6 to 2.5 keV 
(the energy range with calibration mis-match between the 
two XIS detectors) reduces this probability to 0.24 with very 
little change in spectral parameters. This is not enough to 
claim detection of a 4th cyclotron line. The F-stat based test 
performed above, looks at what fraction of the variance of 
the data-set can be explained by the model. This test looks 
at the percentage of data variance that a new model (with 
all its components) can explain as compared to the percent¬ 
age of data variance explained by the old model (with all its 
components). Unaccounted residuals exist in the soft part of 
our fit spectrum below 3 keV . This might account for the 
low reduction in percentage variance accounted for by the 
addition of a cyclotron line at the higher energies (11 keV ) 

. We checked each of the other RXTE observations individ¬ 
ually with this F-test (results listed in Table |T]|. The broad¬ 


band continuum is not completely sampled by the RXTE 
data alone and this could explain the low improvement in 
for the individual RXTE data-sets. To improve the statis¬ 
tics, we did a joint fitting of all the data-sets using both the 
3 CRSF and the 4 CRSF models. In doing this fit, we let 
the line energies for each observation be independently esti¬ 
mated by leaving them untied. This gave an F-statistic value 
of 1.07 and a probability of chance occurrence of 0.15 (or a 
confidence level of 85 %). 

Finally, we tried a more robust numerical evaluation of 
significance from Monte-Carlo simulations. We did this by 
finding the probability of false detection of the 4th CRSF, 
assuming the 3 CRSF model to be true. By simulating spec¬ 
tra using the XSPEC Irt script, we generated a large number 
(7438) of simulated data-sets following the continuum model 
with 3 CRSFs modified by statistical noise. The Irt script 
generates these data-sets by using the fit covariance matrix 
to make a random draw of the fit parameters. The model 
so obtained, is convolved with the response matrices of the 
individual detectors, and statistical (Poisson) noise is added 
to each such simulated data-set. We searched for the pres- 
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ence of the 4th cyclotron line in each of these data-sets by 
trying to fit them with the continuum and 4 cyclotron lines, 
and compared it to a fit with the continuum and 3 cyclotron 
lines. We tabulated the fit value for each such effort. To 
check if the observed 4th CRSF was significant, we compare 
the statistic of each of the simulation runs against the statis¬ 
tic of the observed data. The form of the statistic we used 
was the F-stat (fst) co nstructed as mentioned before (see 
also lSartore et aDl201,'Tl i. We obtained three instances where 
the simulated fst was as high as the observed fst. However, 
in none of these instances, the 4th CRSF fit with centroid 
near 11 keV . This gave a PCI of 3 in 7438, or a significance 
of 3.5cr for the observed CRSF at 11 keV . 

This result, though needs to be treated with care be¬ 
cause of the poor x^ of th® original fit (731/472) v/s 
(688/469). We proceed assuming the 4 CRSF model to be 
the better one, but only after noting that the poor initial fit 
could influence the Monte-Carlo results. 

We then looked at the validity of the 4 CRSF model in 
describing other data-sets. All the individual RXTE spectra 
do not have wide-band coverage and have lower spectral 
resolution than say the combined Suzaku / RXTE data-set. 
This could be a reason for the marginal improvement in 
the x^ statistic. The fact that multiple observations showed 
signs of 4 CRSF features, with an improvement in x^ would 
lead to an increased relevance of this detection. Additionally, 
given that each of these observations gave nearly similar 
centroid energies for all 4 CRSFs (see Fig. El) would also 
increase the significance of our detection. 

The possible reasons for two observations not requiring 
4 CRSFs to describe their spectra are discussed herewith. 
The Swift and INTEGRAL simultaneous observations not 
fitting with the 4 CRSF model could be either because the 
line at 15 keV does not exist in this observation, or because 
the line is too weak to be detected by the JEM-X detec¬ 
tor due to its much lower collection area and poorer signal 
to noise than the RXTE/PCA. This is where we find the 
high signal to noise of the combined Suzaku / RXTE data 
to be critically important. The RXTE standalone observa¬ 
tion which did not show 4 CRSFs, though had a high signal 
to noise all through the expected CRSF energy ranges. We 
take the detection of only 2 CRSFs in this observation to 
be a valid result and discuss the possible reasons for non 
detection in this data-set in the next section. 

A useful indicator for evaluating the strength of the 
absorption lines, and to see if they are physically relevant 
is the line equivalent width. We computed the equivalent 
width as 
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Figure 5. Line Distribution diagram. This figure shows the vari¬ 
ation in energy, FWHM and equivalent widths of all the lines 
we obtained in each of our observations. The central point repre¬ 
sents the line energy and equivalent width with horizontal bars 
depicting the FWHM and vertical bars representing errors on the 
equivalent width. We use black to indicate the 11 keV line set 
and gray to indicate the 15 keV line set. 



pE2 

EW (keV) = / (1 - AE 

J El 

P{E) = 




{E - E,y,y + 


( 1 ) 

( 2 ) 


where P{E) is the line profile that we use, W is FWHM 
of the line, E^yc is its energy and r is its depth. By simple 
error propagation of the variances obtained from the fit co- 
variance matrix, we can get an idea of errors on the estimate. 
We calculated this for all 4 CRSFs to check for their rele¬ 
vance. The results of such calculations are plotted in Fig. El 
We try to make sense of these results in the next section. 


Figure 6. 4 CRSFs and their trends with changing luminosity. 
Top panel shows variation in CRSF energy with luminosity for 
each of the lines in 4 CRSF model. Bottom panel shows variation 
in equivalent width for the same set of lines. These plots are made 
for all data-sets which can be described using 4 CRSFs. 


4 DISCUSSION 

4.1 Summary of results 

In this work we have analysed the spectra from different X- 
ray observatories covering the 2011 outburst of the HMXB 
4U 0115-1-63. These include eight observations spread across 
twenty days including two wide-band observations compris- 
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Figure 7. Trends in the fundamental lines.Top panel shows 
the effect of not using 4 CRSF lines, which leads to the anti¬ 
correlation with luminosity as reported in previous works. Bottom 
panel shows the change in line equivalent width of the fundamen¬ 
tal of each set of the 4 CRSF model. The vertical line is placed at 
the same luminosity level as in Fig[T] The rising phase data-set, 
which does not seem to follow the trend is marked separately. See 
text for details. 

ing of data from multiple satellites (Swift + INTEGRAL 
and Suzaku -\- RXTE). The use of multiple instruments cov¬ 
ering different energy ranges enabled us to properly model 
the broadband continuum, which in turn yielded better con¬ 
straints on the CRSFs. The two major results from our anal¬ 
ysis are: 

(i) Broadband modelling of the continuum: The 

continuum was well constrained with a blackbody compo¬ 
nent modelling the soft X-ray bands and a cut-off power- 
law for the harder bands. The blackbody component was 
found to be statistically significant in both of the wide¬ 
band observations, as well as the standalone RXTE observa¬ 
tions. Excess emission features in the shape of a Gaus sian, 
as reported in previous works e.g. iMiiller et al.l (l2013ll was 
not required for our analysis. Such Gaussian features have 

J ireviously been attributed to possib l e cyc lotron emission 
Becker fc Wolfll l2007l : iFerrigno et al.l l2014l i , however such 
results were inconclusive. 

(ii) Two sets of cyclotron lines: From the spectral 
analysis we find evidence for presence of two independent 
sets of cyclotron lines, each with a harmonic, with funda¬ 
mental energies centred at 11 keV line and ~ 15 keV line. 
Previous works find either one of the two fundamental lines, 
with the 15 keV line found at lower lum inosities, and the 
11 k eV line found at higher luminosities (iTsveankov et ahl 
l2007lj . resulting in the inference of anti-correlation of the 
line energy and luminosity. To check whether the 4 cyclotron 
lines are two independent sets of harmonics we compared the 
variation of equivalent widths of the lines evaluated following 
eq. H We also investigated variations in line energies with 
the outburst luminosity (summarised in Figures O [5] and [T]). 
We took the line equivalent width of the non-detected line 
to be zero. From the figures we note that: 


(a) There is no luminosity dependence of the funda¬ 
mental cyclotron line at 11 keV if we use the 4 CRSF 
model, where both the 11 and ~ 15 keV lines are 
present. The anti-correlation appears if the 3 CRSF model 
is used, as shown in the top panel of Fig. 0 This could 
be an artifact resulting from incorrect spectral modelling 
with only 3 cyclotron lines, instead of 4 CRSFs, as we 
discuss below. 

(b) In the 4 CRSF model, the line equivalent widths 
of the two lines at 15 keV and 11 keV show opposite 
variations with luminosity as seen in the bottom panel 
of Figure [T] 

When only 3 cyclotron lines are used to model the spec¬ 
tra, the two lines at ~ 11 keV and ~ 15 keV are modelled 
by a single CRSF component. Since the equivalent widths of 
these two lines changes with luminosity, the single averaged 
CRSF component is closer to the line with higher equiva¬ 
lent width in the given observation. Simulations have shown 
that often the second harmonic is deeper and more promi - 
nent llArava-Gochez fc Hardin j2000l : [Schonherr et al]|2007lj 
due to photon filling and emission by de-excitation near the 
fundamental CRSF energy. Result s from previous outbursts 
llLi et al.ll2012l : [Boldin et al.ll2013h and our analysis demon¬ 
strate the near constant line energy of the 20 keV line. This 
result too would lead us to expect a near constant line en¬ 
ergy for the fundamental 11 keV line. 

If we split the 4 CRSFs that we obtained into two sets of 
harmonics with one at (11 keV and 20 keV ) and the other 
at (15 keV and 33 keV ), it helps us explain our observations 
as listed below : 

(a) As shown, this explains the reason for observations 
of an anti-correlated fundamental CRSF with luminosity. 

(b) It could be the reason why this source is the only 
one to have shown multiple (upto 5) harmonics of the 
fundamental cyclotron line in its spectrum. If we have two 
such line forming regions with fundamentals at 11 keV 
and 15 keV , and each region showed the presence of 2 
harmonics, then they can easily be confused for multiple 
harmonics from a single 11 keV fundamental CRSF. 

(c) This would explain why the CRSFs we obtain at 
20 keV varies so litle, whereas the one we get at 33 keV 
varies a lot more. Under our hypothesis, the ~ 33 keV 
line in our observations would be a combination of the 
2nd harmonic from the 15 keV set (at ~ 30 keV ) and 
the 3rd harmonic from the 11 keV set (at ~ 33 keV ). As 
seen from Fig. [5] and [O] the ~ 33 keV line shifts to higher 
energies when the ~ 11 keV line becomes stronger. 

(d) Finally, it gives a plausible reason for detection of 
only two CRSFs in the RXTE observation taken during 
the rising phase. If, for some reason, the 11 keV line set is 
either not present or has very weak signatures in the rising 
phase of the outburst, then the observed CRSF energies 
in the rising phase would correspond to the 15 keV line 
set only. The values we get for the observed CRSFs seem 
very close to this. 

4.2 Possible origin of the cyclotron line sets 

In this section we discuss where the two scattering regions 
could be located and what causes the CRSF strengths from 
these regions to vary in such a manner. 
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The first possibility is emission from two different re¬ 
gions at different heights on the same pole. If we assume a 
dipole like magnetic field structure of the Neutron Star (NS) 
the difference of ~ 4 keV between two sets of lines can be 
caused if the difference in height of the emitting regions is 
~ 1.1 km, which is appr oximately the shoc k height above 
the neutron star surface llBecker et al.ll2012h . If the 15 keV 
line set originates in the fan beam at the base of the mound, 
and the 11 keV line set in the pencil beam from the top 
of the shock, then variation in line equivalent widths with 
luminosity can be explained in terms of varying strengths of 
the fan and pencil beam emissions. At higher luminosities, 
the fan beam is expected to dominate, whereas both beams 
will be visible at interm ediate luminosities levels (see Fig. 1 
of I Becker et aLlfioi^ i. However, the effect of such varia¬ 
tions of luminosity on the equivalent width is unclear and 
not adequately addressed by existing theoretical models. 

The second po ssibility is emission from different poles. 
ISasaki et al.l (120121 ') have analysed the pulse profiles of this 
source and found a ~ 60° offset between the position of one 
of the poles and the antipodal position of the other pole. 
This non dipolar field structure could lead to different local 
fields at the scattering regions of each pole, which in turn 
could cause the difference in energies of the lines formed 
in each of them. The variation in accretion rate onto one 
pole as compared to the other can cause the variation in the 
equivalent widths as seen in bottom panel of Fig. 0 

We can get further indicators by looking at the pulse 
resolved energy spectrum and pulse phase-lag spectrum. We 
examined the e nergy dependant pha se lag of the pulses 
as computed by iFerrigno et al.l l|201ll i for the observation 
with simultaneous RXTE/Suzaku data. We have computed 
the phase lags using the data from the higher time reso- 
lution RXTE PC A d etector, in a manner similar to that of 
IFerrigno et al.l (120111 1 . As seen in Figure [H the pulse has two 
distinct peaks. We obtain the phase lags by taking phases 
corresponding to the main peak at a reference energy and 
cross-correlating them against the same set of phase bins at 
other energies. We have performed this for both the main 
and secondary peaks, with phase bins taken as shown in 
bottom panel of Fig. [8l In IFerrigno et all (1201 il l, the au¬ 
thors report that the most negative phase shifts occur at 
energies near the CRSF energies. This has been attributed 
to a change in beam pattern at the CRSF energies, with 
the photons at the CRSF band lying along a pencil beam 
whereas the rest of the emission is dominated by a fan beam. 

In our analysis, we found that main peak’s phase-lag 
spectrum had a similar set of minima at ~ 11 keV , ~ 23 keV 
and ~ 39 keV . However, the secondary peak’s phase-lag 
spectrum had minima at ~ 16 keV and ~ 30 keV . The 
low count rate and increasing errors for energies above 45 
keV made it diflicult for us to find these dips at higher en¬ 
ergies. The phase bins for performing the cross-correlation 
were chosen in such a manner that the pulse profile at the 
reference energy band (t aken as the band betw een 8.17 keV 
to 10.63 keV , similar to IFerrigno et all (I 2 OIII II would have 
a prominent peak like structure in these phase bins. If we re¬ 
fer to the decomposed pulse profiles in Fig. 8 of ISasaki et al.l 
(l2012tl for the decay phase of the outburst, we see that the 
two peaks in our overall pulse profile roughly correspond to 
the emission from the two poles. This lends support to the 
two pole possibility. However, we cannot rule out the fan 


beam / pencil beam as we have not examined the possi¬ 
ble pulse profile and phase lags caused by such a possibil- 
ity. Modelling these “wavy” phase lags, as at tempted before 
(IFerrigno et al.lI 2 OIII : ISchonherr et al.l 120141 1 could help re¬ 
solve between these two possibilities. Pulse phase resolved 
spectra could have given additional indicators. However, the 
large time resolution of the Suzaku XIS data (of 2s) pre¬ 
vented us from getting the spectra at small phase bins of 
the pulsar spinning at 3.6s. 

One way to confirm our hypothesis would be to use a 
single wide-band large area detector to make similar ob¬ 
servations during the next outburst of this source. The 
constraint of using multiple instruments to bridge over 
the energy region having CRSFs and the region having 
the blackbody spectrum makes it very important to have 
these instruments cross-calibrated properly. While the cross¬ 
normalisation constants that we use compare favourably 
with the calibration carried out bv lTsuiimoto et al.l (1201 il l 
for all instruments except Suzaku PIN-XIS, we do note that 
uncertainties in instrument cross-normalisations lead to un¬ 
certainties in the computed flux values and line equivalent 
widths estimated. Having fewer number of instruments to 
cover the range would then be ideal. For example, five to 
six snapshot observations over different luminosities using 
the XMM Newton and NuSTAR telescopes would definitely 
give a higher signal to noise data and lesser uncertainty in 
order to confirm or reject our hypothesis. Another way to 
do this would be the possible construction of a polarization 
spectrum. It is a well known fact that cyclotron resonant 
scattering has highly enhanced cross-sections for incident 
light polarized in the direction parallel to the local magnetic 
field, versus its cross-section for light polarized perpendicu¬ 
lar to it. This enhancement occurs at the res onant scattering 
energ y and drops off at other energies (see iBecker fc Wolfll 
l2007l . and references there-in). So, potentially a measure¬ 
ment of polarization in small energy bands can pin-point 
the energy range over which resonant scattering occurs, and 
effectively delineate the presence of actual CRSFs from those 
resulting due to incorrect spectral modelling. Future obser¬ 
vations by detectors prop osed for measuring the polarization 
in dif ferent energy bands (iPaul et al.ll201ol : Havashida et al.l 
I 2 OI 4 II will definitely help improve our understanding of this 
problem. 


4.3 Conclusion 

In this paper, we demonstrate the utility of having wide¬ 
band high signal to noise data by making use of Suzaku and 
RXTE satellite data. Using this, we point out the presence 
of a blackbody component. We also note the possible indi¬ 
cations of two sets of CRSFs in this source. We note that 
having two sets of lines at 11 keV and 15 keV explains 
the reason for the observed anti-correlation in the funda¬ 
mental CRSF energy with source luminosity. It additionally 
explains the reason for this source being the only known 
accretion powered pulsar to show 5 harmonics of the fun¬ 
damental cyclotron line. Data from large area wide-band 
telescopes like NuSTAR and XMM or from telescopes which 
can give a polarization spectrum can potentially help solve 
this long open problem. 
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Figure 8. Phase lag spectrum and the pulse profile. Bottom panel shows the pulse profile with shaded regions depicting the phase-bins 
over which correlation was computed. Top left panel is the phase-lag spectrum (and the correlation coefficient) for the Main peak (shaded 
dark gray in the pulse profile). Top right panel is the phase-lag spectrum from the secondary peak (shaded light gray in the pulse profile). 
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APPENDIX A: FITTING THE CONTINUUM 
BLACKBODY 

The blackbody component that we use for the continuum 
modelling does not consider effects of compt on scattering of 
this c omponent from the accreting plasma. iFarinelli et al.l 
(l2012l) introduced an XSPEC model for implementing this. 
The results of using this comptonized blackbody are detailed 
in this section. 

There are nine parameters required to describe the 
comptonized blackbody against two required for the sim¬ 
ple blackbody. For getting the model to fit, we had to freeze 
the plasma parameters to the values obtained by previous 
attempts to model this source using such a bulk comptoniza- 
tion of seed photons. 0. This model fits for the plasma elec¬ 
tron temperature (kTe = 1.3 keV ) and optical depth (r 
= 0.41), velocity profile of the accreting plasma (two pa¬ 
rameters, T] = 0.5, P = 0.22), radius of accretion column 
(fo = 0.1) and alb edo percentage (A = 1) from the NS sur- 
face(see Table 1 of IFarinelli et al.]l2012l ). The remaining free 
parameters are the blackbody temperature and the model 
normalization. The results of such a fit are summarized in 
Table EB As seen on comparing these results with Table 
[3l we get blackbody seed temperatures to be similar to the 
model with a simple blackbody. However, the radius of the 

® http://adsabs.harvard.edu/abs/2013A%26A..553A.103F 
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Table Al. Blackbody temperature and radius for wide-band ob¬ 
servations. 


Instrument 

Day (MJD) 

kT ( keV ) 

Radius (kms)^ 

Swift / INTEGRAL 
Suzaku / RXTE 

55736.34 

55750.82 

1 

-^•-^^+0.12 

'^+0.01 

6-67;i|« 

fi 91 —115 


“ calculated from the normalization parameter 

blackbody drops to more reasonable values, although with 
larger errors. This gives a strong indication that a blackbody 
component is indeed required to model the spectrum. 
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